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The supersymmetric Quantum Mechanics approach is applied to embed bound states in the 
energy gaps of periodic potentials. The mechanism to generate periodicity defects in the first Lame 
potential is analyzed. The related bound states are explicitly derived. 
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Let H = —-§£1 + V{x) be a Schrodinger operator with 
V(x) a real-valued and locally bounded function on R. 
The spectrum a(H) of H, as well as its resolvent set 
p{H), will be then complementary subsets of K. The 
involved eigenvalue equation Hip(x) = eip(x), admits 
solutions for arbitrary e € K. The e-values such that 
e = £6 ad{H) lead to square integrable wave- functions 
4>e(x) in L 2 (R), while for e G p{H) one gets eigen- 
functions ip(x) deprived of physical meaning. It is well 
known that the class of Schrodinger's operators H admit- 
ting analytic expressions for the corresponding eigcnfunc- 
tions is narrow. The supersymmetric (susy) Quantum 
Mechanics pj intends to expand the class. In partic- 
ular, the conventional Infold- Hull factorization method 
has been extended to employ not only the physically in- 
terpretable ips but also the nonphysical eigenfunctions 
xp 01 ■ The main idea of this generalized factorization 
is to implement the Darboux-Crum transformation as 
the mechanism to construct multi-parametric families 
{V(Ai, A2, A„, 2;)} of (almost) isospectral exactly solv- 
able potentials 0, 3 • The members of each one of these 
families share common analytical properties, such as the 
number and type of (possible) singularities as well as the 
asymptotic behaviour. Hence, it is intended that the 
procedure does not add new singularities nor modifies 
the asymptotic behaviour of the initial potentials. This 
mechanism is ensured by taking eigenfunctions ip with 
e G p(H) as the key transformation functions Q|. 

The present contribution deals with the generalized 
factorization method applied to periodic potentials. As 
we shall see, the very interesting phenomenon of peri- 
odicity defects in crystal lattices can be produced as a 
consequence of the susy mechanism. We support our 
analysis mainly on the results reported recently in p|. 

Let us consider a Schrodinger operator Hq for which 
the solutions u(x, e) to the eigenvalue equation 

Hou(x,e) = eu{x,e) (1) 

are known for any eel. Notice that equation does 
not require the square integrability of u(x,e). An im- 
portant example is the Hamiltonian with a real-valued, 
continuous and periodic potential 

V {x + t) = V (x), xeR. (2) 

In this case, the translation operator Tu(x, e) = u(x + 
t, e) , maps the two-dimensional subspace H. e of solutions 



of the eigenvalue equation Q into itself. As a conse- 
quence, there is a basis of eigenfunctions {vi(x),V2(x)} 
in Tit for which T admits a 2 x 2 matrix representation [T] 
(the monodromy matrix). The eigenvalue equation for 
[T], evaluated at x = T, determines two eigenvalues t± 
for a given energy e, i.e., t± — t±(e). The corresponding 
eigenfunctions v±(x, e) exist for any e S R and are called 
Block functions. The structure of v± (x, e) depends essen- 
tially on the values of e, namely, if e is such that t± G S 1 , 
then e £ o-(Hq) and the related Bloch functions become 
complex and bounded. The spectrum of Ho is then con- 
tinuous, composed by a series of spectral bands in M.. The 
edges of that bands are given by the e-values such that 
t± = ±1, and the Bloch functions become r-periodic or 
T-antiperiodic. Finally, if e is such that t± £ R/{±1}, 
then e £ p(Ho) and the related Bloch functions are real 
and diverge at either —00 or +00. Due to this mecha- 
nism, p(Hq) splits into a (finite or infinite) sequence of 
forbidden gaps. In summary, ct{Hq) = {e ; t± £ S 1 }, 
while p(H ) = R/a(H ). 

As regards to the supersymmetric treatment of peri- 
odic Hamiltonians, one can consider two different op- 
tions: (i) The Darboux transformations are such that, 
when applied to systems described by J3J), the new po- 
tentials are also r-periodic. (ii) The procedure does not 
respect the global periodicity of the initial potential. In 
the former case, Dunne and Feinberg [(| have shown that, 
sometimes, the nonsingular first order transformed po- 
tential becomes simply a half-period displaced copy of the 
initial one. Such potentials were called self-isospectral. A 
further study (using the band edge eigenfunctions) pro- 
duced the 2-susy self-isospectral potentials |7} ■ The above 
results have been recently generalized to embrace arbi- 
trary J-displaced copies of the initial periodic potential 
[5j . Such a phenomenon has been named the translational 
Darboux invariance (or simply Darboux invariance), and 
it reproduces the self-isospectral case for 5 — r/2. The 
Darboux invariance is achieved by using Bloch solutions 
as transformation functions for e either at the band edges 
or in the resolvent set of the periodic Hamiltonian Hq. 
Concerning the option (ii), we first remark that arbi- 
trary linear combinations of v±(x, e) for e G p(Hq), when 
used as transformation functions u(x, e) can produce reg- 
ular non-periodic modifications in the periodic potential 
@- As we are going to show, moreover, the functions 
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u(x, e) allow to embed discrete energy levels in the gaps 
of Hq and to determine the corresponding square inte- 
grable eigenfunctions 

The generalized first order factorization method |2j 
allows to obtain the solutions </v (x, e) of a trans- 
formed eigenvalue equation Hx(e)(f) e r(x,e) — e'(f> e i(x,e); 

Hi(e) := -51+^1(^1 e), where the new potential Vi(x, e) 
is related with the initial one Vq(x) by means of a Dar- 
boux transformation 

V x {x,e) = V o (x) + 20'(x,e). (3) 

The function j3 here is nothing but the general solution 
of the Riccati equation 

-f3'(x,e)+p 2 (x,e) = V (x)-e. (4) 

At this stage, it is interesting to notice that the trans- 
formation P(x,e) = — -j^lnu(x, e) allows to simplify the 
expressions for the potential (0) and the involved eigen- 
functions <f> t i {x, e) leading to: 

V 1 (x,e) = V (x)-2^lnu(x,e), 



where, by simplicity, we have written u = u(x, e) and 
Vv = ipe>{x). The symbol W(u, ip e >) represents the Wron- 
skian of u and ?/v . 

The analytical behaviour of <fi £ i (x, e) strongly depends 
on the values of e and e'. If e' = E G c(H ) and E > e G 
p(H ), with E G <r(H ) denoting the lowest band edge, 
the eigenfunction (f>E(x,e) can be bounded. A particu- 
lar detail of the method is that the function 4> e (x,e) oc 
l/u{x, e) satisfies identically Hi(e)(f> e (x,e) = e0 e (x, e) 
and becomes orthogonal to all the 4>e{x, e). Moreover, it 
can be chosen in i 2 (M). Thereby a(Hi(e)) = a(H )ll{e}. 
In other words, the Hamiltonian -ffi(e) can be con- 
structed in order to have a bound state embedded in 
the lowest spectral gap of H . Hence, H 1 (e) will be al- 
most isospectral with Ho for a properly selected e. This 
situation allows to consider them as susy partner Hamil- 
tonians. Concerning the analytical properties of V\{x, e), 
it is very common to find that those conditions on e, 
making (j)E{x,e) a square integrable function, lead to po- 
tentials J2J with exactly the same number and nature of 
divergences as the initial one (see for example Q)- In 
particular, as for e < Eq the function u(x, e) can be cho- 
sen nodeless, the potential Vi(x, e) in (JSJ is just a regular 
non-periodic modification of Vo{x). The Figure 1 shows 
one of these first order regular deformations for the sim- 
plest Lame potential 

V (x) = 2msn 2 (x\m). (5) 

Let us remark that the functions u{x, — 0.1), produc- 
ing the potential V\{x,— 0.1) in Figure 1, preserve the 
smoothness of (0 as well as its asymptotic behaviour. 
Moreover, the deformation of V\ in the neighbourhood of 



x = indicates the presence of the bound state energy 
level at e = —0.1 in the first forbidden gap (— 00, Eo). For 
the present case, the related wave-function l/u{x, — 0.1) 
is a smooth function of x, sharply peaked at x = 0. 
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FIG. 1: The Lame potential Q with m — 0.5 (gray) and 
its first order regular deformation (black) with a bound state 
energy level at e = —0.1 < Eo = 0.5 (dashed). 

The procedure can be now iterated by considering e = 
e a as given and V\(x, e a ) as the known exactly solvable 
potential. Thus, one has to look for a new potential 
V2(x,e a ,e) satisfying the second order relationship 

V 2 {x,e a ,e) = Vi(x,e a ) + 2p'(x,e a ,e). (6) 

Therefore, the new function (3{x,e a ,e) must fulfill the 
Riccati equation —(3'{x, e Q , e)-\-0 2 {x, e a , e) = V\{x, e a ) — e, 
where e is arbitrary. The solutions of this equation are 
easily determined by knowing those of for at least two 
different factorization constants e, e a |4j. Then, one can 
use the finite-difference Backlund algorithm It is re- 
markable that (3{x, e a , e) appears to be symmetric under 
the interchange of e a by e and vice versa. This means that 
one can achieve the potential © by finding first V\{x, e a ) 
and then V^rr, e a , e) or, first V\{x, e) and then Vi{x, e, e a ) 
0. More details concerning further iterations and sym- 
metry properties of the Darboux transformations can be 
found in [9| and references quoted therein. Observe that 
potential © and the resulting eigenfunctions 4> e >(x, e a , e) 
can be rewritten as: 

V 2 (x, e a , e) = V (x) - 2^ \nW(u a , u); 

6 ,(x e e) oc w ( u -" u ^.') 
9e'(x, e a , t) oc W (u aiU ) 

where, for the sake of simplicity, we have written u a = 
u{x, e a ). Let us consider the precedent 1-susy step 
as given with e a < Eq G g(Hq), The Hamiltonian 
Hi(e a ) is therefore regular and (f> e '(x, e a , e) becomes an 
eigenfunction of H2(e a ,£) with eigenvalue e'. If e' = 
E G cr(iJ ), e ^ e a and E > e G p{H ), one easily 
finds that 4>E{x,e a ,e) can be bounded. Yet, the func- 
tion (j} e (x,e a ,e) oc l/0 £ (x, e Q ) is also orthogonal to all 
the 4>E(x,e a ,e) and can be chosen in L 2 (IR). Hence 
a(H2(e a ,t)) — <j{Hi(e a )) U {e} and the Hamiltonians 
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Hi(e a ), H2(e a ,e) are almost isospectral. The iteration 
chain now establishes that if 2 (e a , e) is also a susy partner 
of Hq of the second order, i.e., H2(e a , e) and Hq are now 
"sususy partners" . In the present case, the iteration pro- 
cedure leads to nonsingular Hamiltonians at each step, 
hence, the sususy transformation is said to be reducible. 



Figure 2 shows the result of an irreducible second order 
regular transformation of potential (JSJ) for m — 0.5. The 
values of the e's are in the gap (Ei,E[) = (1, 1.5). On 
the other hand, Figure 3 illustrates one of the related 
bound states for e = 1.45. 
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Fig. 2. The initial Lame potential (gray) and its second order 
regular deformation (black) allowing two bound state energy 
levels at e a = 1.1 and e = 1.45. Observe the asymptotic 
translational invariance of the potentials. 



Now, let us consider the case when Eq < e a E p(Hq). 
As the function u(x, e a ) in JSJ has a finite number of ze- 
ros in (0,t), the potential Vi(x,e a ), as well as the func- 
tions (j) e i(x,e a ), become singular. However, there is still 
a chance to construct regular potentials V2(x,e a ,e) and 
functions 4> 6 >(x, e a , e). It is enough to take e a ,e G p(Hq) 
such that Ej < e a ,e < Ey, with Ej,Eji two subsequent 
band edges of H (we have used the notation of 0). 
Hence, for an appropriate choice of the involved parame- 
ters, the Hamiltonian i?2(ea, e) is regular although the in- 
termediary ones Hi(e a ) and Hi(e) are singular. We then 
say that the second order transformation is irreducible. 
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Fig. 3. The probability density for one of the bound states 
associated to the potential of Figure 2 at e = 1.45 



Finally, let us remind that the presence of bound states 
in the continuum is not a novelty in Quantum Mechan- 
ics. For example, Stillinger and Herrick used a modifi- 
cation of the von Neumann- Wigner approach to produce 
local potentials with bound eigenstates embedded in the 
continuum of scattering states |l0j |. Similar results are 
obtained by means of the susy Quantum Mechanics • 
However, the situation when the bound states are embed- 
ded in the energy gaps of a periodic potential has been 
out of the attention for researches in the area. Further 
details on this subject can be found in p|. 

The support of CONACyT project 32086E and of CIN- 
VESTAV project JIRA'2001/17 is acknowledged. 



[1] F. Cooper, A. Khare and U. Sukhatme, "Supersymme- 
try in Quantum Mechanics" , World Scientific, Singapore 
(2000). 

[2] B. Mielnik, J. Math. Phys. 25 (1984) 3387; D.J. 

Fernandez C, Lett. Math. Phys. 8 (1984) 337. 
[3] C.V. Sukumar, J. Phys. A 19 (1986) 2297; A. A. Andri- 

anov, M.V. Ioffe, F. Cannata and J. P. Dedonder, Int. 

J. Mod. Phys. A 10 (1995) 6989; V. M. Eleonsky, V.G. 

Korolev, Phys. Rev. A 55 (19997) 2580; V.G. Bagrov, 

B.F. Samsonov, Phys. Part. Nucl. 28 (1997) 374; B.F. 

Samsonov, Phys. Lett. A 263 (1999) 274. 
[4] D.J. Fernandez C, Int. J. Mod. Phys. A 12 (1997) 171; 

D.J. Fernandez, M.L. Glasser and L.M. Nieto, Phys. Lett. 

A 240 (1998) 15; J.O. Rosas-Ortiz, J. Phys. A 31 (1998) 

L507; J. Phys. A 31 (1998) 10163. 
[5] D.J. Fernandez, B. Mielnik, O. Rosas-Ortiz, B.F. Sam- 
sonov, Phys. Lett. A 294 (2002) 168; D.J. Fernandez, B. 

Mielnik, O. Rosas-Ortiz, B.F. Samsonov, J. Phys. A 35 



(2002) 4279. 

[6] G. Dunne, J. Feinberg, Phys. Rev. D 57 (1998) 1271. 

[7] D.J. Fernandez, J. Negro, L.M. Nieto, Phys. Lett. A 275 
(2000) 338; J. Negro, L.M. Nieto and D.J. Fernandez C, 
Czech. J. Phys. 50 (2000) 1303. 

[8] D.J. Fernandez C.,V. Hussin and B. Mielnik, Phys. Lett. 
A 244 (1998) 309; B. Mielnik, L.M. Nieto and O. Rosas- 
Ortiz, Phys. Lett. A 269 (2000) 70. 

[9] L.M. Nieto, B. Mielnik and O. Rosas-Ortiz, in Pro- 
ceedings of the "VII Encuentro de otofio de Geometria 
y Ffsica", p. 237, M.A. del Olmo and M. Santander 
(Eds), Real Sociedad Matematica Espanola, Madrid, 
Spain (2001). 

[10] F.H. Stillinger and DR. Herrick, Phys. Rev. A 11 (1975) 
446. 

[11] J. Pappademos, U. Sukhatme and A. Pagnamenta, Phys. 
Rev. A 48 (1993) 3525. 



